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We examine the scaling law B oc M a which connects organismal metabolic rate B with organismal 
mass M, where a is commonly held to be 3/4. Since simple dimensional analysis suggests a — 2/3, 
we consider this to be a null hypothesis testable by empirical studies. We re-analyze data sets for 
mammals and birds compiled by Heusner, Bennett and Harvey, Bartels, Hemmingsen, Brody, and 
Kleiber, and find little evidence for rejecting a = 2/3 in favor of a — 3/4. For mammals, we find 
a possible breakdown in scaling for larger masses reflected in a systematic increase in a. We also 
review theoretical justifications of a = 3/4 based on dimensional analysis, nutrient-supply networks, 
and four-dimensional biology. We find that present theories for a = 3/4 require assumptions that 
render them unconvincing for rejecting the null hypothesis that a — 2/3. 



I. INTRODUCTION 



The "3/4-law" of metabolism states that organismal 
basal metabolic rate, B, is related to organismal mass, 
M, via the power law §, |, |, §, g § 



B = cM c 



(1) 



where a is believed to be 3/4. The assumption that a = 
3/4 is relevant in medicine @, |, g], nutrition |lC], |TT|, [l2) , 
and ecology [g, [h| [l4|, [l|, and has been the subject of a 
series of theoretical debates @, ^ [jj], |(], |lj . It has 
been oft quoted that quarter-power scaling is ubiquitous 
in biology M, ^|. Such quarter-law scaling reinforces, 
and is reinforced by, the notion that basal metabolic rate 
scales like B oc M 3 / 4 . 

Nevertheless, the reasons, biological or otherwise, for 
why a = 3/4 have remained elusive and their elucida- 
tion stands as an open theoretical problem. A recent 
surge of interest in the subject, including our own, has 
been inspired by the elegant attempt of West, Brown and 
Enquist p2|] to link nutrient-supply networks to metabol- 
ic scaling. This work suggests that a fundamental under- 
standing of the relationship between basal metabolism 
and body size is within our grasp and that closer inspec- 
tion of both theory and data are duly warranted. 

In this paper we work from the null hypothesis that 
a = 2/3. In a resting state, heat is predominantly lost 
through the surface of a body. One then expects, from 
naive dimensional analysis, that basal metabolism scales 
as surface area which scales as V 2 ^ 3 where V, volume, is 
proportional to M presuming density is constant. This 
scaling of surface area with mass has found strong empir- 
ical support in organismic biology g, 23 Such a 



surface law of metabolism was first expounded in the 
nineteenth century |p5f . Later observations of devia- 
tions from a = 2/3 eventually led to its replacement by 
a ~ 0.72-0.73 which was then supplanted by the sim- 
pler a = 3/4 p|, p6|. The widespread agreement 
that a — 3/4 is due largely to the formative influence of 
Kleiber [g, gj, ^7| and has been accepted and used as a 
general rule for decades 

We re-examine empirical data available for metabol- 
ic rates of homoiothcrms as well as carefully review 
both recent and historical theoretical justifications for 
a = 3/4. Our statistical analysis of data collated by 
Heusner [^8| for 391 species of mammals and by Ben- 
nett and Harvey for 398 species of birds shows that 
over considerable, but not all, ranges of body size, the 
hypothesis a — 2/3 is not rejected by the available data. 
We also review empirical studies by Bartels |50|, Hem- 
mingsen (^3), Brody |2(|, and Kleiber and find the 
data, upon re-examination, to be supportive of our inter- 
pretations. We then examine theoretical attempts to con- 
nect metabolic rate to mass. These include approaches 
based on dimensional analysis [H , |l8| , pi] , |32] , |33| , pi[ , 
four-dimensional biology |3q , p6[ , and nutrient-supply 
networks |22|, |37]]. We find that none of these theories 
convincingly show that a = 3/4, rather than a = 2/3, 
should be expected. 

However, we do not suggest that the 3/4-law should be 
replaced by a 2/3-law of allometric scaling. We instead 
argue for a more general approach to the subject, using 
a = 2/3 as a null hypothesis which should be tested by 
empiricists when considering a ^ 2/3, and acknowledge 
the possibility of deviations from simple scaling. 



II. MEASURING THE METABOLIC 
EXPONENT 
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The history of metabolic scaling may be traced through 
a series of heavily cited empirical papers, some of which 
are composed of very few data points. In order to bet- 
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FIG. 1: Metabolic rate, B (watts), as a function of mass, 
M (kg), for 391 species of mammals. Data taken from 
Heusner [ psf . The straight line represents the best fit for 
the 357 species with mass less than 10 kg where the a — 
0.668 ± 0.019. The upward deviations for species with larger 
mass (see Table may indicate a real biological difference 
but may also be due to the paucity of data. 



ter understand the scaling of metabolic rate, we work 
back in time, calculating a and deviations from uniform 
scaling for data from Heusner [^8) , Bennett and Har- 
vey [£9|, Bartels J|(|, Hemmingsen p3| , Brody f2o| , and 
Kleiber ||. These papers represent some of the most 
influential, widely cited, and often controversial papers in 
the field. Our re-analysis of the data demonstrates that 
a = 2/3 should not be rejected for mammals with mass 
less than approximately 10-35 kg, and a similar analysis 
of metabolic data for birds demonstrates a = 2/3 should 
not be rejected for birds in general. 

We have used the same methods to calculate a and 
its dependence on M in all cases where data is avail- 
able. Slopes and intercepts are determined using stan- 
dard linear regression in log-space taking M to be the 
independent variable. The standard correlation coeffi- 
cient is denoted by r while that obtained using the Spear- 
man rank ordering |38| is written as r s . When data is not 
available we have attempted to classify the data sets in 
terms of the original calculations of a and its dependence 
on M. 
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°"( Q ) 
0.441 
0.093 
0.038 
0.016 
0.010 
0.009 
0.009 
0.008 
0.008 



0.05 
[-0.811,1.719] 
[0.545,1.034] 
[0.578,0.778] 
[0.620,0.704] 
[0.643,0.693] 
[0.651,0.698] 
[0.675,0.720] 
[0.686,0.728] 
[0.689,0.731] 



0.01 
[-1.294,2.202] 
[0.473,1.106] 
[0.550,0.806] 
[0.608,0.716] 
[0.636,0.700] 
[0.645,0.705] 
[0.668,0.727] 
[0.680,0.734] 
[0.684,0.7371 



TABLE I: The exponent a measured for varying ranges of 
mass, M < A/ max (kg), for Heusner's data J28| . For each 
mass range, N is the sample number and the errors, a(a), are 
for two standard deviations. The 95% and 99% confidence 
intervals of a are listed in the last two columns. For small 
mammals (M < 0.01 kg, 17 species) a large error is apparent 
but for increasing M max , « centers around 2/3. A gradual 
upwards drift in a is evident for M max > 10 kg. 



A. Heusner (1991) 

Data on basal metabolic rate for 391 mammalian 
species compiled by Heusner Q is reproduced in Fig- 
ure |l|. Heusner proposed that species could be separated 
into two groups, one of animals whose basal metabolism 
is normally distributed about a regression line and one of 
statistical outliers. Both groups were found by Heusner 
to satisfy a 2/3-law for metabolism. 

The results of simple regression analysis over various 
mass ranges for Heusner's data are shown in Table Q. We 
observe a break in scaling occurring at around M ~ 10 
kg. For those ranges with an upper mass M max < 10 
kg, a — 2/3 appears to be robust. Note that the data 
comprises 179 species of the order rodentia ranging over 
more than three orders of magnitude of mass from 0.007 
kg to 26.4 kg. On separating out these samples, we still 
find a = 0.675 ± 0.025 for the remaining species with 
M < 10 kg and a = 0.681±0.035 for the rodentia species. 

Upon addition of mammals with mass exceeding 10 kg, 
the exponent steadily increases. Given the small number 
of samples of large mammals, one can only speculate on 
the reason for this possible deviation. Primarily, it may 
indicate a real upwards deviation from scaling, with larg- 
er organisms actually having greater metabolic rates than 
predicted by a = 2/3 jl6], |3(], [3{J. Larger organisms are 
reported to scale allometrically in form so such a devia- 
tion may be a result of changes in body shape and hence 
surface are a ffl ^]. Support for this notion comes from 
Economos ]3l|] who finds the relationship between mam- 
malian head-and-body length and mass is better fit by 
two scaling laws rather than one. He identifies 20 kg as 
a breakpoint, which is in accord with our findings here, 
suggesting that geometric scaling holds below this mass 
while allometric quarter-power scaling holds above. The 
upper scaling observed by Economos could also be viewed 



3 



as part of a gradual deviation from geometric scaling. 

The upwards shift of metabolic rates for larger mam- 
mals could otherwise point to problems of measure- 
ment (note the corrections for elephants in Brody's 
data p6[ ), an evolutionary advantage related to larger 
brain sizes jS], |dj , or the lack of competition for ecolog- 
ical niches for large mammals creating a distinction with 
smaller mammals. 



B. Bennett and Harvey (1987) 



Birds show strong support for not rejecting the null 
hypothesis a — 2/3. Figure || shows metabolic data for 
398 distinct bird species taken from Bennett and Har- 
vey n, Hi- We find here that a = °- 666 ± °- 013 

in agreement with Bennett and Harvey's calculations. 
Lasiewski and Dawson |43[ ] similarly found that a = 
0.668 for a smaller set of data. Attempts to reconcile the 
3/4-law with these measurements have centered around 
the division of birds into passerine (perching birds) and 
non-passerine species (non-perching birds). Lasiewski 
and Dawson, for example, found exponents 0.724 and 
0.723 for passerine and non-passerine species respective- 
ly. Though this is not an arbitrary division (core tem- 
peratures are thought to differ by 1-2°C), later work by 
Kendeigh |Q finds exponents ranging from 0.668-0.735 
when passerines and non-passerines are grouped accord- 
ing to different measurement conditions (winter vs. sum- 
mer, etc.). 

Similar distinctions between intra- and inter-species 
scaling have been raised in the study of metabolic scaling 
for mammals where it has been suggested that a = 2/3 
for single species comparisons and a = 3/4 holds across 
differing species fj], §, [jlj. Bennett and Harvey || also 
found that a depends on the level of taxonomic detail one 
is investigating. It remains unclear whether such subdivi- 
sions reflect relevant biological distinctions or underlying 
correlations in the choice of taxonomic levels. 



C. Bartels (1982) 

Bartels |3^] analyzes a set of approximately 85 mam- 
malian species. Although data is not provided in the 
paper, a summary of his results can be found in Table ||. 
Bartels finds a = 0.66 (no error estimate is given) for 
mammals with mass between 2.4 x 10" 3 kg and 3800 kg 
and concludes that the deviation from the expected 3/4 
scaling is due to the variations in metabolic rates of small 
animals. This lends further weight to our conjecture that 
there may be a mass dependence of metabolic rate scal- 
ing. 
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FIG. 2: Metabolic rate, B (watts), as a function of mass, 
M (kg), for 398 species of birds taken from Bennett and Har- 
vey [E9| . The straight line is the result of a regression analysis 
which gives a = 0.666 ± 0.014. 



N 



2.4 x W 5 
2.4 x 10~ 3 
0.26 



3800 ~ 85 0.66 0.99 
0.26 ~ 40 0.42 0.76 
3800 ~ 45 0.76 0.99 



TABLE II: Exponents measured for varying ranges of mass 
(kg), M m ; n < M < M m ax according to Bartels feoj. Here N 
is the sample number and r is the correlation coefficient. 



D. Hemmingsen (1960) 

Hemmingsen's data set |23j for mammals comprises 15 
data points with masses between 0.01 kg and 3500 kg. 
Most of his data is derived from earlier work by Brody. 
He states that the data is well modeled by a power law 
with a = 0.73. To reach this conclusion he does not com- 
pute the power law of best fit, but rather, the "straight 
line . . . was chosen corresponding to [a =]0.73, as estab- 
lished by Kleiber and also by Brody." 

Hemmingsen also finds that a 3/4-law holds for uni- 
cellular organisms. Hemmingsen's work has been cited 
extensively in support of the claim that the 3/4-law is a 
universal biological phenomenon [|[ M, [| ^2). A careful 
re-examination of Hemmingsen's work by Prothero |45f| 
showed that a can range from approximately 0.60 to 0.75 
depending on which unicellular organisms are included in 
the regression. In addition to these questions about seal- 
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0.98 



TABLE III: Results of regression on Brody's data J26|| over 
different mass intervals, A/ m j n < M < Af ma x. An increase in 
a occurs for ranges over larger masses. Here, N is the number 
of data points and r„ is the Spearman correlation coefficient. 



for the data to be 0.738 ± 0.016. Again we consider the 
possibility of a crossover and separate the data into a set 
of 5 species with M < 10 kg and 8 species with M > 10 
kg. For M < 10 kg, a = 0.667 ± 0.043 while for M > 10 
kg, a — 0.754 ± 0.048. These results are again consistent 
with our assertion of a mass-dependent a. Nevertheless, 
it is important to remain mindful of the relative paucity 
of data in these influential works. 



III. FLUCTUATIONS ABOUT SCALING 



ing for unicellular life, Patterson |16| ^7] has theoretically 
shown for aquatic inverterbrates and algae that the scal- 
ing exponent can range from 0.31 to 1.00 depending on 
the mass transfer mechanisms involved. We agree with 
Prothero's conclusions that "a three-quarters power rule 
expressing energy metabolism as a function of size in uni- 
cellular organisms generally is not at all persuasive" E5| . 



E. Brody (1945) 

One of the most influential works on scaling and 
metabolism is that of Brody |26| . Indeed, the scaling law 
for metabolism is sometimes cited as the Brody-Kleiber 
law. Brody compiles a list of metabolic rates for 67 mam- 
mals. The complete data set yields a = 0.73 ±0.01. 
However, on inspection, one makes the surprising obser- 
vation that 32 data points are artificial in that most of 
these are calculated using previously determined empiri- 
cal equations while a few have been corrected to account 
for variations in animal activity. Using the remaining set 
of 35 animals we nevertheless find ev = 0.72±0.02. It is 
also important to note that Brody's research was done 
before the widespread use of electrocardiography and, as 
pointed out by Kinnear and Brown [Q, Brody's data 
may contain overestimations of basal metabolic rate. 

We re-analyze Brody's raw, uncorrected data for mam- 



mals over different mass ranges as shown in Table III 



Again, an increase in a is observed for ranges of larger 
masses. This is consistent with the results from Heusner's 
and Bartel's data which suggest a deviation from perfect 
scaling with increase in mass. Furthermore, it is evident 
that a = 0.72 as calculated by regression on the full data 
set is misleading. We reiterate that we are not suggesting 
that there is any robust scaling law for large masses. The 
results of the regression analysis merely suggest a depen- 
dence of a on the mass ranges being considered and that 
a strict power law may not be appropriate. 



F. Kleiber (1932) 

In his now-famous paper on metabolic rate, Kleiber 
analyzed 13 species of mammals with average mass rang- 
ing from 0.15 to 679 kg We find the scaling exponent 



The next logical step after measuring the metabolic 
exponent and systematic deviations thereof is to consider 
fluctuations about the mean. This is seldom done with 
power law measurements p9| and researchers concerned 
with the predictive power of a scaling law for metabolic 
rate have often pointed to organisms that deviate from 
predictions as being either problematic or different 
|26| , p8| , Hj| . We take the view that fluctuations are to be 
expected and quantified appropriately. 

We thus generalize the relation B = cM a by consider- 
ing P(B \M), the conditional probability density of mea- 
suring a metabolic rate, B, given a mass, M, 



P{B \M) = M- a f{B/M a ) 



(2) 



where the leading factor of M a is for normalization and 

Jo°° f( x ) dx = !• 

Our null hypothesis is that fluctuations are Gaussian 

in logarithmic space, i.e., / is a lognormal distribution. 
Gaussian fluctuations are typically assumed in statistical 
inferences made using regression analysis [ p0[ . Demon- 
strating that / is not inconsistent with a normal distri- 
bution will therefore allow us to use certain hypothesis 
tests in the following section. 

If equation (ph is correct then the sampled data can be 
rescaled accordingly to reconstruct /, the scaling func- 
tion. To do so, one must first determine a. We suggest 
the most appropriate estimate of a corresponds to the 
case when the residuals about the best fit power law are 
uncorrelated with regards to body mass. This is similar 
to techniques used in the analysis of partial residuals ]5lf| 
and we make use of it later. We obtain residuals for the 
range 0.5 < a < 1.5 where the prefactor c of B = cM a 
is determined via least squares. The Spearman correla- 
tion coefficient r s is then determined for the residuals and 
recorded as a function of a. We then take the value of 
a for which r s = as the most likely underlying scaling 
exponent. 

We find r s = when a ~ 0.665 for mammals using 
Heusner's data with M < 10 kg and r s ~ —0.42 when 
a = 3/4. For the entire range of masses in bird data of 
Bennett and Harvey, r s = when a ~ 0.671. 

With these results in hand, we extract / for mam- 
mals and birds, the results for mammals being shown 
in Figure ||. We find the form of / agrees qualitative- 
ly with a lognormal. In order to quantify the quality 
of this agreement we employ the Kolmogorov-Smirnov 
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range 



P 



mammals M < 1 0.153 0.232 0.120 0.307 
mammals M < 10 0.153 0.093 0.120 0.135 
birds all 0.132 0.032 0.115 0.573 

TABLE IV: Results from Kolmogorov-Smirnov tests for the 
underlying distributions of fluctuations around pure scaling 
for both mammals |28| and birds ps| ]. The distribution is 
assumed to be lognormal, i.e., normal in logarithmic coordi- 
nates. The standard deviations a are calculated directly from 
the residuals themselves and determine a level of significance 
p. The a* correspond to p* , the maximum p- value possible. 



test |pC}| , a non-parametric test which gives a signifi- 
cance probability (p-value) for whether or not a sam- 
ple comes from a hypothesized distribution. Not having 
a hypothesis for the value of the standard deviation, we 
take two approaches to deal with this problem. Asserting 
the measured sample standard deviation a to be that of 
the underlying normal distribution, we calculate the cor- 
responding significance probability, p. Alternatively, an 
estimate of a, a*, may be obtained by finding the value 
of er which maximizes p such that p(cr*) = p* . Results for 



both calculations are found in Table IV . All p- values are 
above 0.01, i.e., none show very significant deviations. 
Additionally, the p-value for only the case of the birds 
using a estimated from the data falls below 0.05 indicat- 
ing its departure is significant, but this is balanced by 
the high p-value found by the maximizing procedure. 

Thus, we suggest the data supports the simple hypoth- 
esis of lognormal fluctuations around a scaling law with 
a ~ 2/3. 



IV. HYPOTHESIS TESTS 

We now construct two types of hypothesis tests to 
determine whether or not a = 2/3 or a = 3/4 should 
be rejected by the available data. The first test is the 
standard method of testing the results of a linear regres- 
sion against a presumed slope. The second is a natural 
extension of examining fluctuations about a linear fit as 
per the previous section. By analyzing the correlations of 
the residuals from the best fit line we are able to quanti- 
tatively determine which values of a are compatible with 
the data. In both tests, we reject an hypothesis when 
p < 0.01. 



A. Comparison to a fixed a 

For a given set of N measurements for both mass, Mi, 
and metabolic rate, Bi, we pose the following hypotheses: 




Hq : a = a , 
Hua^ a'. 



(3) 
(4) 



log 10 fi/M- 



FIG. 3: Fluctuations in metabolic rate for mammals with 
M < 10 kg taken from Heusner's data set [Q. The scaling 
function / (see equation (^)) is fitted with a lognormal distri- 
bution. Values of B have been renormalized as B/M 2 ^ 3 and 
partitioned into 20 bins. 



We test the null hypothesis, Ho, in the specific cases a' = 
2/3 and a' = 3/4 for data from KleiberH, Brody ||§, 
Bennett and Harvey [ p9[ , and Heusner pq , over various 
mass ranges. Here, the p-value represents the probability 
that, given two variables linearly related with slope a' 
and subject to Gaussian fluctuations, a data set formed 
with N samples would have a measured slope a differing 
at least by \ot — a'\ from a' [^0). For a null hypothesis 
with a = a', we write the p-value as p a i, e.g., P3/4. 

For mammals with M < 10 kg, the results of the 
hypothesis test are contained in Table |v|. The null 
hypothesis that a — 3/4 is rejected for both Brody and 
Heusner's data and should not be rejected in the case of 
Kleiber. The alternative null hypothesis that a = 2/3 is 
not rejected for both Heusner and Kleiber and rejected 
in the case of Brody. Again, divisions into mass ranges 
are somewhat arbitrary and are chosen to help demon- 
strate the mass-dependence of a. For example, for mam- 
mals with M < 1 kg, Brody's data implies we should not 
reject the hypothesis that a = 2/3. 

Table VI details results for mammals with M > 10 
kg. In the smaller data sets of Kleiber and Brody the 
hypothesis that a = 3/4 is not rejected while for the 
larger data set of Heusner, a — 3/4 is rejected. In all 
cases the hypothesis that a = 2/3 for large mammals is 
rejected. Even though Brody and Kleiber's data sets are 
consistent with an exponent a > 3/4, the relative lack 
of metabolic measurements on large mammals and the 
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TV a aja) P2/3 P3/4 

Kleiber 5 0.667 0.016 0.99 0.088 

Brody 26 0.709 0.020 < 10~ 3 < 10" 3 

Heusner 357 0.668 0.010 0.91 < 10~ 15 

TABLE V: Hypothesis test based on standard comparison 
between slopes that a = 2/3 and a = 3/4 for mammals with 
A/ < 10 kg. Here, a is the measured exponent, a is the stan- 
dard error, and the p- values P2/3 and P3/4 for the hypothesis 
a — 2/3 and a = 3/4 are listed in the last two columns. 



strong rejection by Heusner's larger sample prevents us 
from drawing definitive conclusions about the particular 
value, if any, of a for M > 10 kg. 

N a cr(a) P2/3 P3/4 

Kleiber 8 0.754 0.021 < 10" 4 0.66 

Brody 9 0.760 0.038 < 10~ 3 0.56 

Heusner 34 0.877 0.088 < 10~ 12 < 10~ 7 

TABLE VI: Hypothesis test based on standard comparison 
between slopes that a = 2/3 and a = 3/4 for mammals with 
M > 10 kg. See Table ^ for the definition of all quantities. 

When all mass ranges are considered fo r bo th birds 
and mammals the hypothesis test (see Table |VIl| ) demon- 
strates that both a = 2/3 and a — 3/4 are rejected based 
on the empirical data on mammals, while a — 2/3 is not 
rejected and a = 3/4 is rejected based on the empirical 
data on birds. In summary, we find that a single exponent 
may be appropriate for rough estimates but, from a sta- 
tistical point of view, it appears that no single exponent 
explains the data on metabolic scaling for mammals. 




2 0.6 2/3 0.7 3/4 0.8 0.6 2/3 0.7 3/4 0.8 




0.6 2/3 0.7 3/4 0.8 0.6 2/3 0.7 3/4 0.8 



FIG. 4: Test of the null hypothesis H : r S jAzi,Mi) = 
based on mammalian data from Heusner ] ]28[ ] (see equa- 
tion (|s|)). Shown are plots of p(a) for differing mass ranges. 
In all plots the two dashed horizontal lines correspond to 
p — 0.05 and p = 0.01. The individual plots correspond 
to the following ranges: (a) M < 3.2 kg, (b) M < 10 kg, (c) 
M < 32 kg, (d) all mammals. For all mass ranges considered, 
P2/3 > 0.05 and p 3/4 < 10 -4 . 



N a a{a) P2/3 P3/4 

Kleiber 13 0.738 0.007 < 10" 6 0.11 

Brody 35 0.718 0.011 < 10~ 4 < 10~ 2 

Heusner 391 0.710 0.008 < 10~ 6 < 10~ 5 

Bennett and Harvey 398 0.664 0.007 0.69 < 10~ 15 

TABLE VII: Hypothesis test based on standard comparison 
between slopes that a — 2/3 and a = 3/4 for birds and 
mammals over their entire mass range. See Table [v| for the 
definition of all quantities. 



B. Analysis of residuals 

As per our discussion of fluctuations, a sensitive test 
of a null hypothesis is to check the rank-correlation coef- 
ficient of the residuals. In order to test the hypothesis, 
a = a' , we pose the following hypotheses: 

H : r s>a ,(zi,Mi) = 0, (5) 
Hx-.r^izuMi) £ 0. (6) 



where the Zi are the residuals. The hypothesis Hq means 
that if the residuals for the power law B — cM a are 
uncorrelated with M then a' could be the underlying 
exponent. The alternative hypothesis H\ means that the 
residual correlations are significant and the null hypothe- 
sis should be rejected. The p- values represent the proba- 
bility that the magnitude of the correlation of the residu- 
als, \r s ^ a r(zi, Mi)\, would be at least its value as expected 
for samples taken from randomly generated numbers. 

In this case we have tested the hypothesis for a range 
of exponents, a' = 0.6-0.8, and calculated the signifi- 
cance levels for both mammal and bird data compiled by 
Heusner |2f| and Bennett and Harvey p9| , over differ- 
ent mass ranges. The results of this hypothesis test for 
Heusner's data is contained in Figure || and for Bennett 
and Harvey's data in Figure |^. Both tests show that the 
hypothesis a — 3/4 is rejected while that of a = 2/3 
is not rejected over all mass ranges considered for both 
birds and mammals. This does not mean that a = 2/3 is 
the "real" exponent, but rather that it, unlike a = 3/4, 
is not incompatible with the data. 
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McMahon |l| have tried to explain why quarter-power 
scaling in body lengths and widths should lead to a = 
3/4. Cross-sections of limbs are argued to scale as M 3 / 4 
and therefore the power required to move scales in the 
same way. However, it is not clear why the power output 
of muscles should be the dominant factor in the scaling of 
basal metabolic rate. Furthermore, such quarter-power 
scaling for animal shape is not generally observed [Q. 

Recent debates have focused on deriving a solely from 
dimensional analysis [54j. The problem with all 
attempts to derive metabolic rate from dimensional anal- 
ysis is that different constraints lead to different choic- 
es of contributing scales [Q. Explaining the scaling of 
metabolic rate is therefore displaced to biological ques- 
tions of energetic constraints, mass density, physiological 
time, and diffusion constants across surfaces. 



B. Nutrient Supply Networks 



FIG. 5: Test of the null hypothesis H : r s>a t(zi,Mi) = 
based on bird data from Bennett and Harvey P9J (see Fig- 
ure^ for details). Here, the individual plots correspond to the 
following ranges: (a) M < 0.1 kg, (b) M < 1 kg, (c) M < 10 
kg, (d) all birds. As for the mammal data, P2/3 > 0.05 and 
P3/4 < KT 4 for all mass ranges considered. 



V. THEORIES 

Thus far we have presented empirical evidence that a 
is mass dependent and that the null hypothesis a = 2/3 
should not be rejected for mammals with M < 10 kg and 
all birds in most available data sets. What then of the- 
oretical attempts to derive the 3/4-law of metabolism? 
We show below that many of these arguments, while 
often elegant in conception and based on simple physics 
and geometry, contain sufficient flaws to render them 
unconvincing for the rejection of the simplest theoreti- 
cal hypothesis, a = 2/3. 



A. Dimensional analysis 

Dimensional analysis is a very useful technique when 
there is only one mass, length, and time scale in a giv- 
en problem. However in the case of metabolic scaling in 
biological organisms there has been a long history of the- 
oretical debates over which scales to use when predicting 
the scaling of metabolic rate via dimensional analysis. 

Theories of biological and elastic similarities have been 
used to explain many structural aspects of organisms 
such as the length and width of major limbs |u| B2l {33| . 
Using the principles of elastic similarity, Bonner and 



Interest in Kleiber's law resurged with the suggestion 
by West, Brown and Enquist (WBE) |^2| that nutrient- 
supplying networks might be the ubiquitous limiting fac- 
tor in organismal form. This remains an appealing and 
elegant idea and stands as one of the boldest and most 
significant attempts at discerning the underlying phys- 
ical mechanisms responsible for quarter-power scaling. 
Although previous work had addressed the problem of 
optimal network structure (521 53, 54, pa], theoretical 



relations between optimal networks and the scaling of 
basal metabolic rate had never been considered. 

The basic assumptions of WBE are i) homoiotherms 
have evolved to minimize the rate at which they dis- 
sipate energy; ii) the relevant energy dissipation aris- 
es from transport through nutrient-supply networks; Hi) 
these networks are space- filling; and iv) all homoiotherms 
possess capillaries invariant in size. From these four 
assumptions WBE derive three important conclusions: i) 
nutrient-supply networks are fractal; ii) these networks 
contain area-preserving branching; and Hi) metabolic 
rate scales with a — 3/4. However, as we show below 
and in Appendices [X] and [b|, the arguments used are 
mathematically incorrect and as a consequence none of 
the above conclusions may be derived from the explicit 
assumptions. Nevertheless, we find the model appealing 
and potentially useful in understanding a number of bio- 
logical issues. Thus we detail below where the errors lie 
to illuminate the path of future work. For clarity, we use 
the same notation as WBE. For each level k in the net- 
work hierarchy one has ATj. vessels each with length Ik and 
radius with k = 1 being the aorta and k — N being 
the capillary level. Related important quantities are rife, 
7fc and the ratios of number, length and radius from 
levels k — 1 to k. 

Central to the theory is the connection of these network 
ratios to metabolic rate. WBE find that nk = n, j3k = ft 
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and 7fe = 7 are all constants independent of k and that 



logn 

log7/3 2 



(7) 



This depends in part on an assumption, which we discuss 
below, that B oc N c where N c is the number of capillaries. 
They also conclude that 



ft 



and 



7 = ?i- 1/3 , 



(8) 



which gives a = 3/4 in equation (Q). Whereas we show 
below that these relations do not arise from an opti- 
mization principle, they do have simple interpretations. 
The first relation corresponds to networks being area- 
preserving via N k rf, = Nk-ir^._ v The second rela- 
tion follows from a space-filling criterion that N k l^. = 
Nk-ilk-i- Whether or not space-filling networks satis- 
fy these conditions has been discussed by Turcotte et 
al. who consider the more general case of side- 

branching networks and arrive at an equivalent statement 
of equation (|7) where the network ratios ft and 7 are to 
be determined empirically as functions of n. 

WBE minimize energy dissipation rate by minimizing 
network impedance using a Lagrange multiplier method. 
Two types of impedance are considered: Poiseuille 
flow |57| and, for the case of mammals and birds, a more 
realistic pulsatile flow [ j58| . 

We use the Poiseuille case to demonstrate how frac- 
tality is not proven by the minimization procedure. The 
impedance is given by 



N 

fe=o k 



N 



fc=0 



(9) 



where Z k is the effective impedance of the fcth level. As 
WBE show, the equations are consistent and Z is mini- 
mized when 



Ik = ftk 



-1/3 



(10) 



However, the calculations do not require these ratios to 
be level- independent, and as a consequence, the network 
need not be fractal. Further details may be found in 
Appendix K. To see why this is true, we observe that 
equations (|9) and ( |l0|) give 



Z k 



lkft k A n k l Z k _ l 



l-Z, 



fc-i. 



(11) 



Thus, Z, the quantity being minimized, is invariant as 
long as Ikft^n^ 1 = 1 for each k. This shows that in this 
setting, a network can have n k varying with k and still 
be "efficient." A finding of fractal networks would have 
provided a derivation of Murray's empirical law which 
essentially states that ft — n -1 / 3 for the outer reaches of 
the cardiovascular system |39| . 

Regardless of these issues, the assumption of Poiseuille 
flow leads to an approximate metabolic scaling law with 
a = 1. WBE suggest that modeling pulsatile flow will 



provide the explanation for a 
now takes the form 



3/4. The impedance 



1/2 



rrrl /2 N k ' 



(12) 



where h k is the thickness of the vessel wall. Howev- 
er, as explained in Appendix the equations given by 
the Lagrange multiplier technique are inconsistent. For 
example, the equations give h k = — rfc/5 which means 
negative wall thicknesses for blood vessels when they are 
by definition positive ^o) . If reasonable modifications 
are made to circumvent this issue, then the equations lead 
to a — 6/7 rather than a = 3/4. 

In order to obtain the scaling a = 3/4 one could aban- 
don the minimization calculation and assume a fractal, 
space-filling, area-preserving network where B oc N c . In 
support of such an assumption, there is good empirical 
evidence that blood systems are well approximated by 
fractals |6l], ^32|, |64| . With regards to the assumption 
that B oc N Cl direct measurements for capillary density 
(N c /M oc M a ~ x ) are reported by Hoppeler et al. pq | 
with exponents for the scaling of capillary density across 
species ranging from —0.21 ± 0.04 to —0.07 ± 0.11 for 
various regions of muscle. These numbers are in keep- 
ing with higher exponents for the scaling of N c with M 
in the range 0.75-1.00, but whether or not B oc N c is 
itself an unproven assumption. It is probably more likely 
that the number of capillaries scales with the maximum 
metabolic rate which is thought to scale with an expo- 
nent closer to unity [[36| . At rest not all capillaries diffuse 
oxygen simultaneously and the limiting factor for basal 
metabolic rate might not be N c . 

A simpler and more recent theory based on the idea 
of networks has been proposed by Banavar et al. |37[ . 
Here, networks fill Z?-dimensional hypercubes that have 
L D uniformly distributed transfer sites. The theory is 
applied to both three-dimensional organisms and two- 
dimensional river networks. For organisms, Banavar et 
al. find blood volume scales as Vb oc L^ D+1 \ Since 
Banavar et al. further assume that B oc L D and that 
V b oc M, they conclude that B oc M D/( - D+1 \ Thus, 
when D = 3, this gives a = 3/4. 

However, transfer sites are assumed to be invariant in 
size and hence L D appears to be proportional to volume 
V consequently M. Thus, both the scalings Vb oc M and 
Vb oc M( D+1 ^ D are used, creating an apparent incon- 
sistency. The scaling of the distance between transfer 
sites and the distinction between Euclidean and non- 
Euclidean length scales could possibly be clarified to help 
resolve the dilemma. Note that Vb oc M is supported 
empirically p7fl. 



C. Four-dimensional biology 

Over two decades ago it was suggested by Blum Q 
that a = 3/4 could be understood by appealing to a sur- 
face law of metabolism in a four-dimensional space. In d 
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dimensions, the "area" A of the hypersurface enclosing a 
e?-dimensional hypervolume scales like A oc V~a~ . When 
d = 4, A oc V 3 / 4 , although how this could be reconciled 
with our three-dimensional world was not explained and 
the theory has been refuted elsewhere |pq| . 

Recently, an attempt by West et al. [|36[ has been made 
to refine and generalize their earlier work on metabolic 
scaling |2^| using an optimization procedure to explain 
how an effective fourth dimension could yield a = 3/4. 
The idea put forward is that organisms have evolved 
to maximize the scaling of the effective surface area, a, 
across which resources are exchanged. The area a and 
the biological volume v are shown to satisfy the relation 

aocv^+^/^\ (13) 

where e a and e v are exponents to be determined by opti- 
mization. West et al. then introduce the relationship 
v = al where I is a characteristic length of the organism. 
With the further assumption that d«M, equation fll3| ) 
then becomes 

ffiK M( 2 + £ »)/( 3 + 6 «+ £ '), (14) 

where e; = e v — e a . With the conditions that < ei, e a < 
1, West et al. find that e a — 1 and e; = 0. Equation Jl4| ) 
then yields a oc M 3 ! 4 . Assuming a oc B, this gives a = 
3/4. 

However, this result contradicts the geometric fact that 
transfer area can maximally scale as volume, i.e., a oc v, 
which gives a = 1. Indeed, this result is obtained by 
optimizing equation ( |l3|) instead of equation (|lj) . Doing 
so leads to e a = 1 and e v — 0, assuming < e a , e v < 1, 
which gives a oc M, i.e., a = 1. In order to reconcile this 
with the results of West et al., we note that the bounds 
< e a , e v < 1 are overly restrictive. For example, e\ = 
— 1 corresponds to the relevant length I being invariant 
with respect to M and, in this case, equation (Q) then 
gives the same scaling as (|l^), namely, a oc M. Thus, the 
contradiction is resolved and the optimization procedure 
is seen to yield a = 1 rather than a = 3/4. 



VI. DISCUSSION 

The possibility that there might be a simple law to 
explain the scaling of metabolic rates still captures the 
imagination of many seeking to understand what Kleiber 
called "the fire of life" It is perhaps for this reason 
that so many researchers, theorists and empiricists alike, 
have struggled to deduce explanations for the deviations 
from the simplest expectation that a — 2/3. 

The shift from a — 2/3 to a = 3/4 began with the ear- 
ly work by Kleiber and Brody who found a ~ 0.72-0.73 
in limited data sets ||, ^o| . Afterwards it was work by 
Hemmingsen [ p3[ and a general consensus among practi- 
tioners |21| that simple fractions would be a more con- 
venient standard that led to the widespread acceptance 



of a — 3/4. Subsequently, a = 3/4 has often been tak- 
en as fact despite the absence of a comprehensive theo- 
ry and contradictory evidence from large literature sur- 
veys. Most prominent among these surveys are those by 
Bartels J3^|, Bennett and Harvey p9], and Heusner [ B8[ , 
which suggest that a depends on body size and taxonom- 
ic level. 

We have re-analyzed some of the most influential 
empirical data sets in the study of metabolic rate scal- 
ing. We have constructed a set of hypothesis tests which 
show that in the data sets of Kleiber, Brody, Bennett and 
Harvey, and Heusner, pure 3/4-law scaling is not present. 
For both mammals with M < 10 kg and all birds we 
are unable to reject the null hypothesis a — 2/3. For 
mammals with M > 10 kg, systematic deviations from 
a = 2/3 appear to be present in all of the data sets, 
the roots of which might simply be a consequence of a 
change in body shape for large mammals or might point 
to a greater evolutionary advantage of large mammals. 

We have also reviewed historic and recent attempts to 
justify a = 3/4 theoretically. Many of the early efforts 
to explain the scaling of metabolic rates via dimension- 
al analysis and other crude scaling techniques have been 
dismissed in the past. Although recent attempts to link 
metabolic rates to network structure are noteworthy they 
do not prove the stated conclusions. Nonetheless, we 
believe that research exploring the role of geometric form 
and the dynamics of growth in constraining the behavior 
of networks might lead to important progress in organis- 
mal biology. 

Stated simply, after a systematic review of the avail- 
able empirical data and theoretical arguments, we find no 
compelling evidence of a simple scaling law for metabolic 
rate, and if it were to exist, we also find no compelling 
evidence that the exponent should be a = 3/4. 
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APPENDIX A: NETWORK OPTIMIZATION 
CALCULATION FOR POISEUILLE FLOW 

We follow the conventions of WBE and consider the 
case of Poiseuille flow as a means to derive Murray's law 
and the fractal nature of nutrient supply networks. The 
impedance of the network is 



E 



8^h 



(Al) 



Minimizing the network's impedance with the Lagrange 
constraints of fixed mass and blood volume along with 
the assumption of a space filling network leads to the 
auxiliary function, 



JV 



F m (r k ,l k ,N k ,M) = 



N 



k=0 k 
N 



X *r 2 k l k N k + J2 *kN k ll + X M M. (A2) 



fc=0 



fe=0 



Taking partial derivatives with respect to lj, Tj and JV. 
we respectively have 

dF m 8a* 



and 



dF m 
dr 3 

dF n 



- KrW + \rfN j +3\ j N j l? = 0, (A3) 
-4 • 8/ilj 



+ \2irrjljNj = 0, 



(A4) 



1 + \-kv% + Xjlf = 0. 



Considering first equation (A4), we obtain 

16ju 
7r 2 rf JV? ' 

Since this holds for all j then 

16/x 



1 = 



ttM^JV* 16/x ' 3 v 



(A5) 



(A6) 



(A7) 



where (3j = rj/r,-_i and 



Nj/Nj-x which demon- 



strates that 



(A8) 



giving us Murray's law [f59| 



After rearranging equation (A3) we obtain 
8a* 



A; 



Xnrj 
31/ 



16/x 



(A9) 



where we have used the form for A obtained in equa- 
tion (A6). Note that derivatives with respect to Nj, 



equation (A5), yield the same expression for Xj given 
above: 



Ai = - 



-1 • 8n 
^fNflf 
8fi 



Mrrj 

*/ 
16/i 



8/i 



J 3 



(A10) 



The three equations (A3), (A4) and QAq ) are therefore 
consistent but redundant. The redundancy can be seen 
to lie in the fact that the auxiliary function F m in equa- 
tion ( A2) can be written in terms of only two variables for 

and £fc = r k /l k . Equation (|A2|) 



each level k: £ k 
thus becomes 



N, 



F m (Zk,(k,M) 



N 

E 

fc=0 
N 



8/i 



AE 7r ^^+E A ^ fc + AMM - ( An ) 



k=0 



k=0 



Thus one is only able to obtain information such as ratios 
of variables rather than exact values for network param- 
eters. 

The scaling of length ratios are explicitly determined 
by WBE's space- filling assumption 



N k q = a 



(A12) 



Thus, even without implementing the minimization pro- 
cedure the space-filling assumption implies 



Ik 



-1/3 



(A13) 



where Tk — h/h -i is the le ngth r atio. Finally, the equa- 



tions (A8), (A9) (or (A10)) and (|A13|) combine to give 



2 ; 2 
l'fc-l 



U-l 



r 4 N 2 l 2 

4 -2 -2 



1/3 )- 4 % 2 (% 1/3 )- 2 



= 1 



(A14) 



so we have Xk = Xq for all k. 

The calculations are seen to be consistent and yield 
Murray's law (59) . Variations with respect to M are more 
subtle since JV = N(M) and provide higher order correc- 
tions. However, one of WBE's crucial results, nk = n, 
i.e. that the network is fractal, has not been reproduced. 
One way to see this is to consider the impedance as 
impedances in series: 



N 

E 

fc=0 



Z k where Z k 



8filk 



(A15) 



11 



Using equations (A8) and (A13) we have that 



Z k 
Zk-i 



,-1/3 



1. 

(A16) 



In other words, the same impedance appears at each lev- 
el. So 



(N + 1)Z N ~ NZ_ 



N 



N 



AT 8 ^ lc 

N — -. oc — 

nr^Nc N, 



(A17) 



since r c and l c are assumed to be independent of mass and 
N c is the number of capillaries. This is true regardless of 
whether or not the structure is fractal. The network has 
to possess branching ratios that collectively maintain the 
same impeda nce fr om level to level (i.e., 7fc//3^rifc = 1 as 
per equation (A16)) but there is no requirement that the 
individual ratios 7^, (3 k and n k be independent of level. 
Moreover, without the result that the network is fractal, 
this minimization procedure no longer yields the power 
law scaling of metabolic rate (see equation ((?])). 



Note that the extra variable of wall thickness, h k , has 
been included in the second term to make it a measure 
of the total volume taken up by the blood system. The 
variable hk must appear in the constraints if the mini- 
mization is to make any sense and the blood volume is 
the only reasonable choice — the blood volume constraint 
becomes a network volume constraint. 

On considering variations of equation ( |B2| ) with 
respect to Tj and hj we obtain 



dF„ 



5 (Ehjp) 



1/2 



dr 



and 



2 \/2TTry 2 Nj 



^— + X2n( rj + h^ljNj = 0, (B3) 



(B4) 

Since the second term of these equations are the same 
we then have an equality between the first terms which 
simplifies to show that 



APPENDIX B: NETWORK OPTIMIZATION 
CALCULATION FOR PULSATILE FLOW 

In the case of Poiseuille flow, WBE find a network 
structure where area preservation is not satisfied {(3k 7^ 
n k 1 ^ 2 ) and, effectively, a = 1 (if n k = n is assumed). 
The intended fix is to properly treat pulsatile flow of 
mammalian blood circulation systems. By doing so we 



should obtain (3k 



-1/2 



and 7fe 



-1/3 



Together 



with the assumption n k = n, this leads to the conclusion, 
N c oc Af 3 / 4 , and assuming B oc N c , it would imply a 3/4- 
law of metabolic scaling. 

The calculation relies on the results of Womersley's 
work on pulsatile flow ]5^, [jfj. Womersley's calculations 
lead to a modification of the Poiseuille impedance. For 
large tubes one has 



pep 



(Bl) 



where c Q = (Eh/2pr) 1 ^ 2 is the Korteweg-Moens velocity, 
E is Young's modulus, h is the thickness of the vessel 
wall, p is blood density and r is, as before, the vessel 
radius |5^, |6(|. This impedance appears to be per unit 
length but it has the correct dimensions showing that 
for a flow with pulsatile forcing in an elastic tube, the 
impedance is independent of the tube length. 

Womersley's impedance suggests a new auxiliary func- 
tion: 

N (Fh n) 1 / 2 
F w (r k ,h k , lk,N k , M) = J2 ^_5/2 



^0 V2nrl /2 N k 

N 

A + h k ) 2 l k N k + AfcJV fc lg + X M M. (B2) 



k=0 



k=0 



h J =--r r (B5) 

This suggests that r k is the distance to the outer wall of 
blood vessels. We are then measuring the blood volume 
as before and we should have had (r k — hk) instead of 
(rk + hk) in the auxiliary function. However it is apparent 
from Womersley's work [^8], Q that r is the radius as 
measured from the center to the inner wall of a blood 
vessel rather than the outer wall. There appears to be 
no reasonable and simple way of including the hk into a 
constraint function and we have an ill-posed problem. 

Nevertheless, we may proceed with the calculation by 
adding an extra assumption that hk = ao r fe where ciq > 0. 



The modified version of equation (B3) now gives A as 



(a £p) 1/2 



V2vr 2 (l 



a ) 2 rlN 2 l 



(B6) 



j ■) 



Since the right hand side is independent of j we must 
therefore have 



(B7) 



and given the space-filling constraint, jj = n,j 1 ^ 3 , we 
obtain 



1 o4 2 -1/3 a4 5/3 
1 = PjUjUj = p-Tlj 



3 '"3 



(B8) 



which gives a relationship between the radius and number 
ratios that is not area-preserving: 



(3j = n] 



5/12 



(B9) 



A further complication here is that the equations 
obtained by setting dF w /dlj = and dF w /dNj = are 
not consistent. 
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As in the case of Poiseuille flow, rik = n is not deriv- 
able. Assuming that n.& = n and using equation ([?]) we 
find that the metabolic exponent should be 

Inn Inn 
a = -J^W = - mn -i/a n -io/i 2 - 6 / 7 ' ( B1 °) 

as opposed to the stated 3/4. 

1/2 

Note that if we had found (3k — nj then the 3/4 
law would have been deduced (again, assuming rik = 
n). Another observation here is that if the Womersley 



impedance is taken together with 0k — n k then we 
find that the minimum total impedance is obtained irre- 
spective of the ratios n^ being equal or not. So, in the cas- 
es of Poiseuille and pulsatile flow a fractal network is not 
necessary for energy dissipation to be minimized. Addi- 
tionally, in the case of a pulsatile flow network, a — 3/4 
cannot be derived from the optimization problem as stat- 
ed. It may instead be derived by assuming an area pre- 
serving, space- filling, fractal network where B oc N c . 
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